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1. INTRODUCTION 

^— i ! 

^jr)' One motivation for studying the production of charmonium in high energy collisions 

is that these particles are among the simplest probes for heavy-ion collisions. In recent 
years, there have been significant developments in heavy quarkonium production, both in 

Q\ [ theory and in experiment. These developments may have important implications for the 
use of charmonium as a probe in heavy-ion collisions. A summary of these developments 

y—i • is presented below. Those interested in more details are referred to a recent review article 
by Braaten, Fleming, and Yuan |]J. 

m ; 2. CHARMONIUM AS A PROBE 

oo : 

The simplest probes for hard processes in heavy-ion collisions are leptons and photons. 
One can argue that the next simplest are the heavy quarkonium states, bottomonium 
' and charmonium. The bottomonium states are too heavy to be produced in abundance 
in heavy-ion collisions, but charmonium states are produced reasonably copiously. The 
Ph. states with the cleanest experimental signatures are the J PC = 1 states J/ip and ip 1 ', 
which decay into e + e~ and fi + fi~. The J ++ states Xcj, J — 1,2, also have reasonably 
good experimental signatures through their radiative decays into ip + 7. 

The reason that these charmonium states are such simple probes is that they are es- 
sentially 2-body systems. The ip, for example, is to a good approximation a bound 
state of a cc pair in a color-singlet state with angular- momentum quantum numbers 3 Si. 
The color wavefunction is (RR + BB + GG) / a/3, the space wavefunction has the form 
V'(r) — R{f)/V^, where R(r) is the radial wavefunction, and the possible spin wave- 
functions are ft, (Tl + It)/v2, and ||. This cc state is the dominant component of the 
wavefunction, but there are also higher Fock state components. There is, for example, a 
small ccg component. Since the overall color state of the ccg must be color- singlet, the 
cc pair must be in a color-octet state. The 8 independent color-octet states are RG, GB, 
BR, RB, BG, GR, (RR - GG)/y/2, and (GG - BB)/y/2. For most observables, the 
contribution from the ccg component of the wavefunction is small. However, for some 
observables, the contribution from the cc component of the wavefunction is suppressed 
and components in which the cc pair are in a color-octet state can be important. 

There are three important features of charmonium that make it simpler than light 
hadrons. For each of these simplifications, there is a theoretical tool that can be applied 
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to charmonium that is not available in the case of light hadrons. These tools give us a 
handle on the effects of gluons in each of the three most important ranges of wavelength. 
These three ranges of wavelengths are much less than, comparable to, and much greater 
than r, where r is the typical separation of the cc pair in charmonium. The first simplifying 
feature of charmonium is that the mass of the charm quark (m c ~ 1.5 GeV) is much larger 
than the momentum scale Aqcd associated with nonperturbative effects in QCD. This 
allows us to use perturbation theory in the running coupling constant a s (m c ) to calculate 
the effects of gluons whose wavelengths are of order l/m c or larger. The second simplifying 
feature is that the c and c are nonrelativistic: the typical value of the relative velocity 
v of the c and c is about \, so v 2 ps |. This allows us to use a relativistic expansion 
in powers of v 2 to organize the effects of gluons whose wavelengths are comparable to r, 
which is on the order of l/(m c v). The third simplifying feature is that the geometrical 
size of charmonium states is significantly smaller than that of light hadrons: the typical 
separation r of the c and c is only about 0.2 fm. This enables us to use the multipole 
expansion in powers of rj A to simplify the interactions of gluons whose wavelengths A are 
much larger than r. By exploiting these three theoretical tools (perturbation theory, the 
relativistic expansion, and the multipole expansion), we hope to understand charmonium 
physics in sufficient detail to make it a useful probe of heavy-ion collisions. 

3. CHARMONIUM PRODUCTION 

We now turn to the problem of charmonium production, where we will focus on the 
production of the J/ip in hadron collisions. We can think of the production of this state 
as proceeding through two steps. The first step is the creation of a cc pair with small 
relative momentum. If their relative momentum is much larger than m c v, they will be 
unlikely to bind and will instead fly apart and form D mesons. The creation of the cc pair 
must involve momenta on the order of m c or larger, because the momenta of the incoming 
particles provide the rest energy of the c and c. The second step in the production of the 
ip is the binding of the cc pair. This step necessarily involves momenta of order m c v or 
smaller, because gluons whose wavelengths are comparable to the size of the bound state 
play a large role in the binding. 

We first consider the creation of the cc pair. One way to produce them is through 
collisions between partons from the colliding hadrons. The simplest such processes are 
gg — > cc and qq — > cc. Parton collision processes such as these always involve virtual 
particles that are off their mass-shells by amounts of order m c . The cross section for 
creating the cc pair can therefore be calculated as a perturbation series in a s {m c ). One 
can also show that the cc pair from such a process is produced with separation of order 
l/m c or smaller, which is much smaller than the size of a quarkonium state. Thus, 
when the cc pair is produced, it is essentially pointlike on the scale of the charmonium 
wavefunction. 

There are also production mechanisms that produce cc pairs that are not pointlike. For 
example, a gluon from one hadron can undergo a quantum fluctuation into a virtual cc 
pair. This pair can make a transition into a real cc pair by exchanging soft gluons with the 
other hadron. The separation of the resulting cc pair is comparable to the wavelength of 
the exchanged gluon. Such processes are particularly important in forward or diffractive 



3 



production. Fortunately, one can minimize the complications from such processes by 
going to a kinematic region where they are suppressed, such as transverse momentum 
much greater than Aqqd- We will therefore restrict our attention from this point on to 
parton collision processes that produce pointlike cc pairs. 

We now consider the binding of the cc pair to form the ip. Assuming that the cc pair is 
produced through parton collisions, all we need to know is the probability for a pointlike 
cc pair to bind to form a ip. The inclusive differential cross section can be expressed in a 
form in which short-distance and long-distance contributions have been factored: 



The sum over n includes all possible color and angular momentum states of the cc pair. 
The short-distance factor is the cross section da for creating a cc pair in the state n, 
which can be calculated as a perturbation expansion in a s (m c ). The long-distance factor 
(0%) gives the probability for a pointlike cc pair in the state n to bind to form the ip. 
The matrix elements (0%) contain all the nonperturbative information that is required 
to calculate the inclusive cross section. 

Until the last few years, most of the predictions for charmonium production in high 
energy collisions were carried out using the color- singlet model ||. This model is a very 
simple prescription for the probability for a pointlike cc pair to form a charmonium meson. 
The factor (0%) is assumed to be nonzero only if the state n is a color singlet with the 
same angular-momentum quantum numbers as the dominant Fock state of the meson, 
which is 3 Si in the case of the ip. The probability for a pointlike cc pair in the color- 
singlet 3 5*i state to form the ip is the same as the probability for the ip to consist of a 
pointlike cc pair. It is proportional to |i?(0)| 2 , where -R(O) is the radial wavefunction at the 
origin. The color-singlet model is very predictive. It gives predictions for ip production in 
all high energy processes in terms of the single nonperturbative factor |i?(0)| 2 . Moreover, 
the factor |i?(0)| 2 can be measured experimentally by the decay rate for ip —>■ e + e~ . Thus 
the predictions of the color-singlet model are absolutely normalized. 

We now consider the production of prompt charmonium in pp collisions. "Prompt" 
means that the charmonium is produced by QCD interactions rather than by weak decays 
of hadrons containing bottom quarks. The cross section for prompt ip production at 
Fermilab's Tevatron has been measured by the CDF collaboration for transverse momenta 
in the range 5 GeV < Pt < 20 GeV . The background from decays of bottom hadrons 
were removed using a silicon vertex detector. The contribution from the radiative decays 
Xcj ^7 was measured and subtracted. The data is shown in Figure 1. The dotted line 
gives the prediction of the color-singlet model at large pt- (The dashed line gives another 
contribution in the color-singlet model that becomes important at lower p?-) The result 
is rather surprising. The measured cross section was found to be about a factor of 30 
larger than predicted by the color-singlet model. 

This enormous discrepancy between theory and experiment is very surprising, because 
the data extends out to large pt where the theoretical analysis is particularly clean. 
There are factorization theorems of perturbative QCD that guarantee that the inclusive 
differential cross section for the production of a hadron at large pt must be dominated 
by fragmentation, the formation of that hadron in the jet initiated by a parton with large 
transverse momentum. In the case of charmonium, large transverse momentum means 
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Figure 1. CDF data on the differential cross section for prompt ^'s that do not come from 
x's as a function of px- The dashed and dotted lines are two contributions predicted by 
the color-singlet model, while the solid line is the prediction of the color-octet mechanism 
with the normalization adjusted to fit the data. 
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Pt 3> rn c and the most important partons are gluons. The inclusive differential cross 
section can therefore be written in the form 

da(^(P)+X) = f dzda(g(P/z)+X)D g ^(z), (2) 

J 

where da is the differential cross section for producing a gluon with momentum P/z and 
D g -Hi>(z) is the gluon fragmentation function, which gives the probability that the jet 
initiated by the gluon will include a ip carrying a fraction z of the gluon momentum. All 
of the dependence on pt resides in the factor da, which can be calculated using pertur- 
bation theory in a s {pr)- At leading order, it comes from the parton process gg — > gg, 
which has a cross section of order ocKpt)- All effects from momentum scales of order 
m c or smaller are factored into D g ^{z). In particular, the nonperturbative effects as- 
sociated with the binding of the cc pair are all included in the fragmentation function. 
The color-singlet model gives a definite prediction for the fragmentation function in (0). 
The leading contribution comes from the parton process g — > ccgg. The resulting frag- 
mentation function is proportional to al(m c )\R(0)\ 2 jvr? c . The dotted curve in Figure 1 
is the prediction obtained using this fragmentation function. Gluon fragmentation gives 
the dominant contribution in the color-singlet model for px greater than about 5 GeV. 
However, the predicted cross section falls about a factor of 30 below the CDF data. 

How can we make sense of this enormous discrepancy between theory and experiment? 
The factorization theorems summarized by (0) are on a very firm foundation. The cross 
section must certainly have this form at the largest values of px that are available. If the 
prediction fails to agree with the data, the problem must lie in our assumption about the 
fragmentation function D g _,^(z). The color-singlet model must drastically underestimate 
the probability for a gluon to fragment into charmonium. 

4. NRQCD FACTORIZATION APPROACH 

How can the probability for the gluon to fragment into charmonium be so much larger 
than predicted by the color-singlet model? There is a simple and natural explanation that 
is based on a new approach to charmonium production developed recently by Bodwin, 
Braaten, and Lepage This approach is called the NRQCD factorization formalism, 
because it makes use of an effective field theory called nonrelativistic QCD. The basic 
idea is to exploit the relativistic expansion for the long-distance factors (0$) in the fac- 
torization formula ([J). The factor (0$) gives the probability for a pointlike cc pair in 
the color and angular-momentum state labelled by n to bind to form a ip. It can be 
expressed as a well-defined NRQCD matrix elements that scales in a definite way with v, 
the typical relative velocity of the c and c. The leading matrix element in v is proportional 
to |i?(0)| 2 and scales like t> 3 . The corresponding term in the factorization formula ([!]) is 
the cross section of the color-singlet model. All other terms have matrix elements that 
are suppressed by powers of v 2 . Thus, if the parton cross sections da((cc) n + X) were 
all comparable, the color-singlet model would indeed give the dominant term in the cross 
section. 

However, the parton cross sections da((cc) n + X) can vary widely in magnitude. They 
can be suppressed not only by powers of a s {m c ) or a s (pr), but also by powers of kinemat- 
ical factors like m c /pT- If the color-singlet model term is suppressed by such a factor, then 
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other terms may be important in spite of their suppression by powers of v 2 . After |i?(0)| 2 , 
the next most important matrix elements are (Og( 3 Si)), (Og( l S Q )), and (Of ( 3 P )), all 
of which are suppressed by v 4 . They give the probabilities for forming a if) from pointlike 
color-octet cc pairs in spin-triplet S-wave, spin-singlet S-wave, and spin-triplet P-wave 
states, respectively. These matrix elements are nonperturbative quantities that cannot 
be related in any simple way to the cc wavefunction. We have no effective means of cal- 
culating them from first principles. These matrix elements must therefore be treated as 
phenomenological parameters to be determined from experimental data. 

Now consider the fragmentation function D g ^(z) that describes the formation of a if) 
from a high-p^ gluon. The color-singlet model term in the fragmentation function has a 
short-distance factor of order a 3 (m c ) and a long-distance factor that scales like v 3 . All 
other terms have long-distance factors that are suppressed by powers of v 2 , but there are 
some for which the short- distance factor is lower order in a s . In particular, the term that 
corresponds to the parton process g — > cc has a short-distance factor that is only of order 
a s . The matrix element in this term is (Of ( 3 Si)), which gives the probability of forming 
a if) from a pointlike cc pair in a color-octet 3 S\ state. The fragmentation function at 
leading order in a s is 

IW*) = ^p^(l-^)(^( 3 ^))- (3) 

This term in the fragmentation function is of order a s v 7 , compared to a 3 v 3 for the color- 
singlet term. It is not immediately obvious which of these two terms is more important, 
but experience tells us that suppression by a power of a s tends to be more effective than 
suppression by a power of v 2 . Thus it is not unreasonable for this color-octet term to 
dominate. 

If the color-octet term is sufficiently large, it could provide an explanation for the CDF 
data on prompt charmonium production JJ. The color-octet term in the cross section at 
large pt is obtained by inserting the fragmentation function (^) into (^). The resulting 
prediction for the differential cross section as a function of px is shown as a solid line in 
Figure 1. Its shape fits the data quite well. The matrix element (Of ( 3 Si)) can be adjusted 
so that the curve also agrees with the data in normalization. The value of the matrix 
element that is required is (Of ( 3 5i)) = 0.014 GeV 3 . Now for this explanation of the data 
to be viable, the matrix element must satisfy an important consistency check. It must 
be small enough to be consistent with suppression by v 4 relative to the corresponding 
color-singlet factor |i?(0)| 2 . Translating both quantities into the probability density for 
a point-like cc pair to bind to form a if), we get 21/fm 3 for the color-singlet 3 Si state 
and 1.8/fm 3 for the color-octet 3 Si state. The latter value is perfectly consistent with 
suppression by a factor of v A . 

We have shown that a color-octet term in the gluon fragmentation function provides 
a plausible explanation for the CDF data on prompt if> production. To make this expla- 
nation convincing, we have to show that this color-octet production mechanism explains 
other aspects of charmonium production. There are two particularly dramatic predictions 
that should be tested experimentally in the near future. The first prediction is that at the 
Tevatron, ip n s at large pr should be almost completely transversely polarized |J. The if>' 
essentially inherits the transverse polarization of the fragmenting gluon. There is an effort 
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underway at CDF to measure this polarization. The second prediction is that prompt 
^>'s should be observable in Z° decay at LEP 0. The color-singlet model predicts a rate 
for prompt %j) production that is too small to be observed. Using the value of (Of ( 3 5'i)) 
obtained by fitting the CDF data, we find that the color-octet mechanism increases the 
production rate by almost an order of magnitude. The predicted branching fraction for 
Z° — > ip+X is 1.4 x 10~ 4 , which is large enough that it should be observable. Preliminary 
measurements of prompt ip production from some of the LEP detectors seem to indicate 
that the branching fraction does indeed have the magnitude predicted by the color-octet 
production mechanism. 



5. CONCLUSIONS 



If we understand charmonum production in simple high energy collisions, it should be 
possible to predict the production rate in heavy-ion collisions. A new framework has been 
developed for describing charmonium production in those kinematic regions where it is 
dominated by parton collisions. The cross section is factored into parton cross sections 
that can be computed using perturbation theory and well-defined matrix elements that 
give the probability for forming the bound state. The matrix elements with the greatest 
phenomenological importance include |i?(0)| 2 , (0 8 ( 3 Si)), (Os^Sq)), and (C 8 ( 3 -Po)) f° r ^ 
and ip' and |i?'(0)| 2 and (Os^Si)) for Xcj- It should be possible to determine these factors 
with reasonable precision by carrying out a thorough analysis of all the available data on 
charmonium production, including fixed target data from pN, ttN and eiV collisions as well 
as data from the high energy pp, e + e~ , and ep colliders. The resulting matrix elements can 
then be used as inputs for predictions of charmonium production in heavy-ion collisions. 

In conclusion, new data from high energy colliders has been driving progress in heavy 
quarkonium physics. There are powerful theoretical tools that can be used to study 
heavy quarkonium. The developments described above have exploited two of these tools: 
perturbation theory and the relativistic expansion. The multipole expansion is another 
powerful tool that has important implications for heavy quarkonium production ||. By 
exploiting these tools, it should be possible to understand heavy quarkonium in sufficient 
detail to allow it to be used as a probe in heavy-ion collisions. 
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